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Abstract 

■I We propose an improved method for hadron-coUider mass determination of new states that decay to a 

^ massive, long-lived state like the LSP in the MSSM. We focus on pair-produced new states which undergo 

three-body decay to a pair of visible particles and the new invisible long-lived state. Our approach is 
— . to construct a kinematic quantity which enforces all known physical constraints on the system. The 

distribution of this quantity calculated for the observed events has an endpoint that determines the mass 
of the new states. However we find it much more efficient to determine the masses by fitting to the 
I— I entire distribution and not just the end point. We consider the application of the method at the LHC for 

various models and demonstrate that the method can determine the masses within about 6 GeV using 
only 250 events. This implies the method is viable even for relatively rare processes at the LHC such as 
^ ^ neutralino pair production. 

' — ' 1 Introduction 

m 
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At hadron colliders the determination of the masses of new particles associated with missing momentum 
signals is very challenging due to the fact that the kinematics of the event cannot be completely reconstructed. 
Hadron colliders collide partons within each hadron, and each parton involved in the collision carries an 
On unknown fraction of the hadron's momentum. Therefore, the center-of-mass (COM) energy and the frame of 

reference of the parton collision are unknown for each event . The problem is further aggravated because one 
CN does not expect any of the new short-lived particle states to travel far enough to create tracks in the detector. 

In extensions of the Standard Model such as supersymmetry (SUSY) or Universal Extra Dimensions (UED) 
there is often a massive, stable, neutral particle that will leave the detector unnoticed, leading to missing 
momentum associated with the production and decay of the new particles required in such extensions. 

For these reasons there has been much work developing techniques to determine the mass of the new 
particles at hadron colliders such as the LHC. Significant information comes from the endpoints of kinematic 
^ invariant distributions. This is illustrated in the simple case that a short-lived state Y undergoes a three-body 

^ decay to a lepton pair plus the escaping neutral particle N (the LSP in supersymmetry), Y ^ + l~ + N. 

For this decay the invariant mass mfi = {ki~ + ki+)'^ has a maximum value equal to the mass difference 

(My - Mn)^ 

maxm;; — My — M]\[. (1) 

New states appear as bumps in the mu distribution where one can read off the mass difference from the 
upper edge of the bump[^ If Y undergoes a two-body decay to an on-shell intermediate state Y ^ X + l^ ^ 
N + 1^ + 1^ , then the shape of the mu distribution will be more like a right triangle with a vertical drop, and 
the maximum mu is given by mf^ — {My — Mjr){Mjr — M^)/M^. These techniques have been extensively 
used to study SUSY in the context of the LHC (see Ref [3] for a study of several models) . Using events with 
four leptons in the final states and missing energy, Ref [T shows how such edges can form a Dalitz-like plot to 
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^ Loop corrections play a role in shifting this endpoint slightly. For a detailed study see Ref [l] . The shape is determined by 
the degree of interference with the slepton, see [2| for examples. 
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determine information about the mass spectra of new states. In short, such edges can accurately determine 
relations between the masses of the unknown particles but not the mass Mjq. The task of determining the 
complete mass spectra is therefore dependent on determining A/^r. 

Much of the work in determining the il/^r in a hadron collider focuses on a cascade of decays. The idea is 
to use events that contain many final states so that one can find enough edges of invariant mass distributions 
to invert the relationships and solve for the masses or SUSY model parameters. For example Bachacou, 
HinchlifFe, Paige [5] use a sequence >g + X2~^^ +^^ + '7^' +^^ + 'Z + Xi involving four new states in 
the event. One can form four invariant mass distributions from these final states, and one has four unknown 
masses. This set of constraints sometimes has multiple solutions. Fitting the shapes of the distributions can 
lift this degeneracy as was shown in Miller, Osland, Raklev (MOR)[B] and Lester [7]. 

Is there a way to find Mjv if there are only three new states involved in the event? Cheng, Gunion, Han, 
Marandella, McElrath (CGHMM) |8] study pair-produced states, Y, that decay via an on-shell intermediate 
state X. An example scenario would be pair-produced X2S where each branch decays via ^2 ^ + ^ 
l~ + l'^ + X° or its conjugate. Their events consist of four leptons and missing energy. They analyze 
each event's kinematics for compatibility with on-shell condition for the assumed topology. To make their 
approach robust against background and finite resolution error, they form distributions and use the shape 
to determine the unknown masses. 

Finally, what can one determine from events which involve only two new states? Cho Choi Kim and Park 
(CCKP)[ni [TU] show how to use the Cambridge transverse mass variable Mt2 of Lester and Summers [TT] 
to find My (which in their case was the gluino mass) assuming a three-body decay to xl and q, q. Their 
example uses about 40000 events where gluinos are pair produced and decay to four jets and missing energy. 
The Mt2 variable is a function x which is an assumed mass of Mpf. One plots the maximum Mt2(x) over 
the 40000 events as a function of x- A kink appears in the function at the correct Mj^ and My □ Using 
this approach, CCKP find My and Mj^ to about ±2 GeV for the case where M^/M_ « 1.3 where 

M+ = My + Mn M_ = My - Mn. (2) 

In this paper we will concentrate on the latter possibility involving the production of only two new 
states. Our particular concern is to use the available information as effectively as possible to reduce the 
number of events needed to make an accurate determination of My and M^. The main new ingredient 
of the method proposed is that it does not rely solely on the events close to the kinematic boundary but 
makes use of all the events. Our method constrains the unobserved energy and momentum such that all 
the kinematical constraints of the process are satisfied including the mass difference, eqQ, which can be 
accurately measured from the II spectrum. This increases the information that events far from the kinematic 
boundary can provide about My and significantly reduces the number of events needed to obtain a good 
measurement of the overall mass scale. Although we develop the method for the case that Y decays via 
a three-body decay to an on-shell final states Y N + l'^ + , its generalization to other processes is 
straightforwarcPj 

In Section [2] we introduce the M2C distribution whose endpoint gives My, and whose distribution can 
be fitted away from the endpoint to determine Aly and AI^ before one has enough events to saturate 
the endpoint. Section [3] estimates the performance for a few SUSY models where we include approximate 
detector resolution effects and where we expect backgrounds to be minimal. Finally we conclude and discuss 
directions for further research. Appendix A discusses the relationship between our distribution and the 
kink in Mt2{x) of CCKP and how this relationship can be used to find M2C in a computationally efficient 
manner. Appendix B provides details of our simulations. 

^For a recent study on situations which lead to kinks using the transverse mass see ref |12l 

^We note that the on-shell intermediate case studied by CGHMM is also improved by including the relationship measured by 
the edge in the II distribution on each event's analysis. The Y decay channel with an on-shell intermediate state X has an edge 
in the II invariant mass distribution which provides a good determination of the relationship maxm^j = {My — ^^x^^^x ~ 
M'^)/Mx- This relationship forms a surface in ,Mx ,My space that only intersects the allowed points of CGHMM's fig 3 
near the actual masses. 
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Figure 1: We consider events with the new state Y is pair produced and in which each Y decays through a 
three-body decay to a massive state N invisible to the detector and visible particles 1, 2, 3, and 4. 



2 An improved distribution from which to determine My 

We consider the event topology shown in fig [l] The new state Y is pair produced. Each branch undergoes 
a three-body decay to the state N with 4-momentum p (q) and two visible particles 1-1-2 (3-1-4) with 
4-niomentum a (/?). The invariant mass mi2 (TO34) of the particles 1-1-2 (3-1-4) will have an upper edge 
from which one can well-determine M^. Other visible particles not involved can be grouped into V with 
4-momentum k. In the analysis presented here, we assume fc = and check that it remains valid for k < 20 
GeV. 

We adapt the concept from Mt2 of minimizing the transverse mass over the unknown momenta to allow 
for the incorporation of all the available information about the masses. To do this we form a new variable M2C 
which we define as the minimum mass of the second to lightest new state in the event My constrained to be 
compatible with the observed 4-momenta of Y's visible decay products with the observed missing transverse 
energy, with the four-momenta of Y and N being on shell, and with the constraint that Af_ — My — Mn is 
given by the value determined by the end point of the mi2 distribution. The minimization is performed over 
the ten relevant unknown parameters which may be taken as the 4-momenta p and q of the states N, and the 
lab- frame collision energy Pq and longitudinal momenta Pz ■ We neglect any contributions from unobserved 
initial state radiation (ISR). Thus we have 

M|c= min (p+a)2 (3) 
subject to the 7 constraints 

{p + af^{q + (3f, (4) 

(5) 
(6) 

(7) 

Although one can implement the minimization numerically or by using Lagrange multipliers, wc find the 
most computationally efficient approach is to modify the Mt2 analytic solution from Lester and Barr [13 . 
Details regarding implementing M2C and the relation of M2C to Mt2 and the approach of CCKP arc in 
Appendix A. 

Errors in the determined masses propagated from the error in the mass difference in the limit of fc = 
are given by 




where SM_ is the error in the determination of the mass difference M_ . To isolate this source of error from 
those introduced by low statistics, we assume we know the correct M_, and one should consider the error 
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Figure 2: The distribution of 30000 events in 5 GeV bins with perfect resolution and no background. The 
three curves represent My = 200 GeV (dot-dashed), My = 150 GeV (dotted) and My = 100 GeV (sohd) 
each with = 50 GeV. Each distribution cuts off at the correct My. 



described in eq(|8| as a separate uncertainty from that reported in our initial performance estimates in the 
next section. 

Because the true p, q, P^, Pz are in the domain over which we are minimizing, M2C will always satisfy 
M2C < My. The equality is reached for events with either mi2 or 77134 smaller than M_, withp^/po = az/oLo, 
and qz/qo = I3z/ Po, and with the transverse components of a parallel to the transverse components of /3. 

The events that approximately saturate the bound have the added benefit that they are approximately 
reconstructed (p and q are known). If Y is produced near the end of a longer cascade decay, then this 
reconstruction allows one to determine the masses of all the parent states in the event. The reconstruction 
of several such events may also aid in spin correlation studies. 

In order to determine the distribution of M2C for the process shown in fig [l] we computed it for a 
set of events generated using the theoretical cross section and assuming perfect detector resolution and no 
background. Details of the simulation are in Appendix B. Figure |2] shows the resulting distribution for three 
cases: My = 200 GeV, My = 150 GeV and My = 100 GeV each with M_ = 50 GeV. Each distribution 
was built from 30000 events. Note that the minimum My for an event is A/_. The endpoint in the three 
examples is cleax, and one is able to distinguish between different My for a given The shape of the 

distribution exhibits only modest model dependency as described in Appendix B. 

One can also see that as M_|_/M_ becomes large, the My determination will be hindered by the small 
statistics available near the endpoint or backgrounds. To alleviate this, one should instead fit to the entire 
distribution. However it is clear that events away from the endpoint also contain information about the 
masses. For this reason we propose to fit the entire distribution of M2C and compare it to the 'ideal' 
distribution that corresponds to a given value of the masses. As we shall discuss this allows the determination 
of My with a significant reduction in the number of events needed. This is the most important new aspect 
of the method proposed here. 
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3 Application of the method : SUSY model examples 



To illustrate the power of the fit to the full M2C distribution, we now turn to an initial estimate of one's 
ability to measure My in a few specific supersymmetry scenarios. Our purpose here is to show that fitting 
the M2C distribution can determine My and Mm with very few events. We include detector resolution 
eff'ects but neglect backgrounds but assume fc = in the simulation. We calculate M2C for the case where 
the the analytic Mt2 solution of Barr and Lester can be used to speed up the calculations as described in 
Appendix A. Details on our calculations and simplifying assumptions can be found in Appendix B. A more 
complete detailed study will follow in a subsequent publication. 

Although fitting the M2C distribution could equally well be applied to the gluino mass studied in CCKP, 
we explore its applications to pair-produced X2- We select SUSY models where X2 decays via a three-body 
decay to l'^ + + Xi- The four momenta a = pi+ + pi- for the leptons in the top branch, and the four 
momenta (3 = pi+ + Pi- for the leptons in the bottom branch. 

The production and decay cross section estimates in this section are calculated using MadGraph/MadEvent 
[H] and using SUSY mass spectra inputs from SuSpect [TS]. The distributions in this section still neglect 
background, but scale the a and j3 four vectors by a scalar normally distributed about 1 with the width of 

^-^4L= + ^ + 0.007 (9) 
ao v/"o(GeV) ao(GeV) 

to simulate the typical LHC detector lepton energy resolution [Tni [T7] . The missing transverse momentum is 
assumed to be whatever is missing to conserve the transverse momentum after the smearing of the leptons 
momenta. We do not account for the greater uncertainty in missing momentum from hadrons or from muons 
which do not deposit all their energy in the calorimeter and whose energy resolution is therefore correlated 
to the missing momentum. Including such effects requires a more detailed detector simulation and is beyond 
the scope of this Letter. These finite resolution effects are simulated in the determination of the ideal 
distribution and in the small sample of events that is fit to the ideal distribution to determine My and M^. 
We do not expect expanded energy resolutions to greatly affect the results because the resolution effects are 
included in both the simulated events and in the creation of the ideal curves which are then fit to the low 
statistics events to estimate the mass. 

We consider models where the three-body decay channel for X2 ^i^^ dominate. These models must satisfy 
rrij^o — m^o < Mz and must have all slepton masses greater than the m^o . The models considered are shown 
in Table [f| The Min-Content model assumes that there are no other SUSY particles accessible at the LHC 
other than ^2 and x° and we place rri-^o and m^o at the boundary of the PDG Live exclusion limit [18]. 
SPS 6, PI, and 7 are models taken from references [H], j2U], and [2T] respectively. Each has the X2 decay 
channel to leptons via a three-body decay kinematically accessible. We will only show simulation results 
for the masses in model PI and SPS 6 because they have the extreme values of M+/M_ with which the 
performance scales. The Min-Content model and the 7 model are included to demonstrate the range of the 
masses and production cross sections that one might expect. 

Bisset, Kersting, Li, Moortgat, Moretti, and Xie (BKLMMX) (3] have studied the 4 lepton -t- missing 
energy standard model background for the LHC. They included contributions from jets misidentified as 
leptons and estimated about 190 background events at a £ = 300 fb~^ which is equivalent to 0.6 fb. Their 
background study made no reference to the invariant mass squared of the four leptons, so one only expects 
a fraction of these to have both lepton pairs to have invariant masses less than A/_ . Their analysis shows 
the largest source of backgrounds will most likely be other supersymmetric states decaying to four leptons. 
Again, one expects only a fraction of these to have both lepton pairs with invariant masses within the range 
of interest. The background study of BKLMMX is consistent with a study geared towards a 500 GeV e+ e~ 
linear collider in ref [22 which predicts 0.4 fb for the standard model contribution to 4 leptons and missing 
energy. The neutralino decay to r leptons also provide a background because the r decay to a light leptons 
I = e, fi (Tr^wjT w 0.34) cannot be distinguished from prompt leptons. The neutrinos associated with 
these light leptons will be new sources of missing energy and will therefore be a background to our analysis. 
The di-r events will only form a background when both opposite sign same flavor rs decay to the same flavor 
of fight lepton which one expects about 6% of the time. 

Table [2] breaks down the LHC production cross section for pair producing two X2 each of these models. 
In the branching ratio to leptons, we only consider e and /i states as the r will decay into a jet and a 
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Model 


Mm Content (ref [18J) 


SPS 6 (ret 1191) 


PI (ref [20 ) 


7 ( ref [2T]) 


Definition 


Xi and X2 
are the only 
LHC accessible 
SUSY States 
with smallest 
allowed masses. 


Non Universal 
Gaugino Masses 
TUo = 150 GeV 
mi/2 = 300 GeV 
tan/3 = 10 
sign(/i) = + 

Ao = 

Ml = 480 GeV 

M2 = A/3 = 300 GeV 


mSUGRA 
TJio = 350 GeV 
mi/2 = 180 GeV 
tan ^ = 20 
sign(^) = + 


Non-Universal 
Higgs Model 
mo = 330 GeV 
mi/2 = 240 GeV 
tan = 20 
sign(/i) = + 
^0 = 

Hi = -(242 GeV) 2 
Hi = +(373GeV)2 


m^o 


46 GeV 


189 GeV 


69 GeV 


95 GeV 




62.4 GeV 


219 GeV 


133 GeV 


178 GeV 




6.6 


13.6 


3.2 


3.3 



Table 1: Models with x'i decaying via a three-body decay to leptons. We only show simulation results for the 
masses in model PI and SPS 6 because they have the extreme values of Af+ /M_ with which the performance 
scales. 









Events with 


Model 


c^xS xl Direct 
cr^o ^0 Via 5 or g 


^^X2^q+q+Xi 


4 leptons +Et missing 
+ possible extra jets 
C = 300 fb^i 


Min Content 


2130 fb 
N/A 


0.067 
0.69 


2893 


SPS 6 


9.3 fb 
626 fb 


0.18 
0.05 


6366 


PI 


35 fb 
12343 fb 


0.025 
0.66 


2310 


7 


17 fb 
4141 fb 


0.043 
0.64 


2347 



Table 2: The approximate breakdown of signal events. 
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Model PI with 250 Events 
Fit Gives M;^,o = 133 +/- 6 GeV 
True M;^.,o=69 GeV, M;^20=133 GeV 




120 125 130 135 140 145 

My o of Theoretical Distribution (GeV) 



Figure 3: fit of 250 events from model PI of Ref [20] to the theoretical distributions calculated for different 
M^o values but fixed M^o - M^o. The fit gives M^g = 133 ± 6 GeV. 



Model SPS 6 with 3000 Events 
Fit Gives My^o= 221 +/- 20 GeV 
True M.,0 = 189 GeV, M.,o = 219 GeV 




180 200 220 240 260 280 

My o of Theoretical Distribution (GeV) 



Figure 4: fit of 3000 events from model SPS 6 of Ref [19] to the theoretical distributions calculated for 
different M^g values but fixed M^g - M^o. The fit gives M^g = 221 ± 20 GeV. 



7 



neutrino introducing more missing energy. Direct pair production of X2 ^ rather modest cross section, 
but production via a gluino or squark has a considerably larger cross section but will be accompanied by 
additional QCD jets. One does expect to be able to distinguish QCD jets from r jets [23^. 

We now estimate how well one may be able to measure m^o and m^o in these models. Figures [Sl and 
[4] show a of the AI2C distribution from the observed small set of events to 'ideal' theoretical Af2C 

distributions parameterized by m^o. The 'ideal' theoretical distributions are calculated for the observed 
value of M_ using different choices for m^o. A second-order interpolation is then fit to these points to 
estimate the value for rtij^o . The 1 a uncertainty for m^o is taken to be the points where the increases 
from its minimum by one. 

The difficulty of the mass determination from the distribution grows with the ratio Af_|_/M_. Figures 
[3] and [4] show the two extremes among the cases we consider. For the model PI M_|_/M_ — 3.2, and for 
model 7 = 3.3. Therefore these two models can have the mass of rriy^o and rriy^o determined with 

approximately equal accuracy with equal number of signal events. Figure [3] shows that one may be able 
to achieve ±6 GeV resolution after about 30 fb~^. Model SPS 6 shown in fig [4] represents a much harder 
case because M+/M_ = 13.6. In this scenario one can only achieve ±20 GeV resolution with 3000 events 
corresponding to approximately 150 fb^^. In addition to these uncertainties, one needs to also consider the 
error propagated from (5M_ in eqQ. 

4 Summary and Conclusions 

We have proposed a method to extract the masses of new pair-produced states based on a kinematic variable, 
M2C1 which incorporates all the known kinematic constraints on the observed process and whose endpoint 
determines the new particle masses. However the method does not rely solely on the endpoint but uses the 
full data set, comparing the observed distribution for M2C with the ideal distribution that corresponds to a 
given mass. As a result the number of events needed to determine the masses is very significantly reduced 
so that the method may be employed at the LHC event for processes with electroweak production cross 
sections. 

We have performed an initial feasibility study of the method for several supersymmetric models. This 
includes the effect of detector resolution but not backgrounds, cuts and combinatoric complications but was 
modeled with an assumption that fc = 0. We demonstrated that for some of the models studied we are able 
to determine the masses to within 6 GeV from only 250 events. This efficiency is encouraging although a 
study including more of the real-world complications is needed to augment this initial study. 

The method we advocate here can be readily extended to other processes. By incorporating all the 
known kinematical constraints, the information away from kinematical end-points can, with some mild 
process dependent information, be used to reduce the number of events needed to get mass measurements. 
We shall illustrate this for other cases elsewhere [in preparation]. 
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Mj2(x) for Four Different Events 
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Figure 5: The Mt2{x) curves for four events with = 50 GeV and My = 100 GeV. Only the events 
whose curves starts off at Mt2(0) > intersect the straight line given by Mt2{x) ~ X — The Mt2 
at the intersection is M2C for that event. 



Appendix A : Using Mt2 to Find M2C 

The variable M7'2, which was introduced in by Lester and Summers is equivalent to 

Mt2{x)= min (p+a)2 

subject to the 7 constraints 

2 2 

p = q 

{Po,0,0, Pz) ^ P + q + a + P + k 

2 2 
P =X ■ 



(10) 
(11) 

(12) 
(13) 
(14) 



As is suggested in the simplified example of |24j , the minimization over Po and Pz is equivalent to assuming 
p and a have equal rapidity and q and /3 have equal rapidity. Implementing this eq( 10 1 reduces to the 
traditional definition of the Cambridge transverse mass. 

By comparing Mt2{x) ^ defined above to AI2C defined in eq([3|, one can see that they are very similar 



with the exception that the constraint eqQ is replaced by the constraint eq(14). M2C can be found by 
scanning Mt2{x) for the x value that such that the constraint in eqQ is also satisfied. 

One can see the M2C and AIt2 relationship visually. Each event provides a curve Mx2 ix) 1 fig [S] shows 
curves for four events with Afjv — 50 GeV and My = 100 GeV. For all events Mt2{x) is a continuous and 
monotonically increasing function of x. As CCKP point out, at large x and at fc = the maximum MT2ix) 
approaches x + so one knows the slope of Mt2{x) for all events will be everywhere less than or equal to 
one. Furthermore if Mt2{x = 0) > M_, as is true for two of the four events depicted in fig|5] then, barring 
an asymptote, there is a solution to Mt2{x) = X + At this point Mt2{x) = min My | constraints = -^20- 
Equivalently 



M2C = Mt2 at X where Mt2{x) 
= A/_ otherwise. 



X + M_ if Mt2{x = 0) > M_ 



(15) 
(16) 



At fc = the maximum x of such an intersection occurs for x — which is why the endpoint of AI2C 
occurs at the correct My and why this corresponds to the kink of CCKP. Because Barr and Lester have an 
analytic solution to Mt2 in ref in the case A; = 0, this is computationally very efficient as a definition. 
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Appendix B: Numerical Simulation Details 



Numerical simulation of "ideal" events 

In order to determine the distribution of M2C for the processes shown in fig [T] it is necessary to generate a 
large sample of "ideal" events corresponding to the physical process shown in the figure. For simplicity in 
the numerical simulations included in this note we always assume fc = and decay via an off-shell Z-boson 
as this is what could be calculated quickly and captures the essential elements to provide an initial estimate 
of our approach's utility. 

Even under these assumptions, one might expect that the shape of distribution depends sensitively on the 
parton distribution and many aspects of the differential cross section and differential decay rates. Surprisingly 
this is not the case; the shape of the distribution depends sensitively only on two properties: 

(i) the shape of the TO12 (or equivalently 77134) distributions. In the examples studied here for illustration 
we calculate the mi2 distribution assuming it is generated by a particular supersymmetric extension of the 
Standard Model, but in practice one should use the measured distribution which is accessible to accurate 
determination. The particular shape of TO12 does not greatly affect the ability to determine the mass of N 
and Y so long as one can still find the endpoint to determine My — Afjv and use the observed mu distribution 
to model the shape of the M2C distribution. 

(ii) the angular dependence of the A'^'s momenta in the rest frame of Y. In the preliminary analysis 
presented here we assume that in the rest frame of X2 1 Xi 's momentum is distributed uniformly over the 47r 
steradian directions. While this assumption is not universally true it applies in many cases and hence is a 
good starting point for analyzing the efficacy of the method. 

Under what conditions is the uniform distribution true? Note that the ^2'^ spin is the only property 
of X2 that can break the rotational symmetry of the decay products. For Xg's spin to affect the angular 
distribution there must be a correlation of the spin with the momentum which requires a parity violating 
coupling. Consider first the Z contribution. Since one is integrating over the lepton momenta, the parity 
violating term in the cross section coming from the lepton-Z vertex vanishes and a non-zero correlation 
requires that the parity violating coupling be associated with the neutralino vertex. The Z-boson neutralino 
vertex vanishes as the Z interaction is proportional to X27^7^Xi^m '-"^ X27^Xi-2^ai depending on the relative 
sign of m^o and m^o eigenvalues. However if the decay has a significant contribution from an intermediate 
slepton there are parity violating couplings and there will be spin correlations. In this case there will be 
angular correlations but it is straightforward to modify the method to take account of correlations. We hope 
to study this in another publicatiorj^ 

Even in the case that the slepton contribution is significant the correlations may still be absent. Because 
we are worried about a distribution, the spin correlation is only of concern to our assumption if a mechanism 
aligns the spin's of the X2S in the two branches. Table [2] shows that most of the X2 '^^^ expects follow from 
decay chains involving a squark, which being a scalar should make uncorrelated the spin of the X2 
branches. One would then average over the spin states of X2 ^^'^ recover the uniform angular distribution 
of x°'s momentum in Xj's rest frame. 

Once one has fixed the dependencies (i) and (ii) above, the shape of the distribution is essentially 
independent of the remaining parameters. To illustrate this result we show in fig [6] two cases: 

(1) The case that the collision energy and frame of reference and angle of the produced Y with respect to 
the beam axis are distributed according to the calculated cross section for the process considered in Section 
[3] in which X2 decays via Z exchange to the three-body state l'^ + + Xii convoluted with realistic parton 
distribution functions. 

(2) The case that the angle of the produced Y with respect to the beam axis is arbitrarily fixed at 6* = 0.2 
radians, the azimuthal angle </> fixed at radians, and the total 4-momentum of the colliding particles 
arbitrarily set to P = (500, 0, 0, 0) GeV. 

The left plot of fig [6] shows the two distributions intentionally shifted by 0.001 to allow one to barely 
distinguish the two curves. On the right side of fig [6] we show the difference of the two distributions with 
the 2 a error bars within which one expects 95% of the bins to overlap if the distributions are identical. In 
addition to tests with fc = 0, we also tested that fc ^ 20 GeV does not change the shape of the distribution 
to within our numerical uncertainties for any of our results. In a test case where we constructed events 

^ Studying and exploiting the neutralino spin correlations is discussed further in Refs [251 1261 127) . 
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Figure 6: Demonstrates the distribution is independent of the COM energy, angle with which the pair is 
produced with respect to the beam axis, and the frame of reference. 



with My = 150 GeV and Mn = 100 GeV, with uniformly distributed with between 2 of 20 GeV, 
with |A:/fco| = 0.98, and with uniform angular distribution, we found the M2C distribution agreed with the 
distribution shown in figure [6] within the expected error bars after 10000 events. Scaling this down to the 
masses studied in PI we trust these results remain unaffected for fc < 20 GeV. Introduction of cuts on jets 
and missing traverse energy will probably introduce some dependence on the COM energy of the collision 
that is absent in this ideal case. 

Given this structure detailed in (i) and (ii) above we calculate the "ideal" distributions for M2C assuming 
that k — and that in the rest frame of Y there is an equal likelihood of N going in any of the Air steradian 
directions. The observable invariant is determined according to the differential decay probability of X2 
to e~ and Xi through a Z-boson mediated three-body decay. Analytic expressions for cross sections were 
obtained from the Mathematica output options in CompHEP 

Inclusion of backgrounds will change the shape. Backgrounds that one can anticipate or measure, like 
di-rs or leptons from other neutralino decays observed with different edges can be modeled and included 
in the ideal shapes used to perform the mass parameter estimation. A more complete study is beyond the 
scope of this letter and will follow in a subsequent publication. 



Least squares fit 

In order to determine My it is necessary to quantify the comparison between the A'' observed events and the 
"ideal" events. To do this we define a distribution by computing the number of events, Cj, in a given 
range, j, (bin j) of M2c- Assuming a Poisson distribution, we assign an uncertainty, aj, to each bin j given 
by 

a]^^-{Nf{M2c„MY) + C,). (17) 

Here the normalized distribution of ideal events is f{M2c, My), and the second term has been added to 
ensure that the contribution of bins with very few events, where Poisson statistics does not appljj^ have a 
reasonable weighting. Then is given by 

.-(M.).i:( ^'""T'"""'-' )'^ (18, 

bin j ^ ^ ' 

The minimum y^iMy^ is our estimate of My . The amount My changes for an increase of by one gives 
our Icr uncertainty, bMy, for My [29]. As justification for this we calculate ten different seed random 
numbers to generate ten distinct groups of 250 events. We check that the My estimates for the ten sets are 
distributed with about 2/3 within ^My of the true My as one would expect for 1 a error bars. One might 

^By this we mean that A'' f{M2Cj, Afy) has a large percent error when used as a predictor of the number of counts Cj when 
N f{M2Cj,MY) is less than about 5. 
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worry that with our definition of x^, the value of per degree of freedom is less than one. However this 
is an artifact of the fact that the bins with very few or zero events are not adequately described by Poisson 
statistics and if we remove them we do get a reasonable per degree of freedom. The determination of My 
using this reduced set gives similar results. 
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